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proach can reproduce effectively the full next-to-leading order numerical 
results derived using non-relativistic QCD (NRQCD) factorization for¬ 
malism. We demonstrate that the next-to-leading power contributions 
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for the rate as well as polarization of J/p production at current collider 
energies. 
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1 Introduction 


Since the discovery of the J/V’, heavy quarkonia, with their clearly separated multi¬ 
ple momentum scales, have been serving as ideal systems to test our understanding 
of QCD bound states and their hadronization processes. Unfortunately, a theoret¬ 
ically and phenomenologically satisfying framework is still not achieved for heavy 
quarkonium production, to explain their yields and the polarizations. The problem is 
more acute with the recently discovered XY Z-mesons, since the production of these 
exotic mesons requires good understanding of the production of conventional heavy 
quarkonia BE]. 

So far the most phenomenologically successful model for heavy quarkonium pro¬ 
duction is based on the non-relativistic QCD (NRQCD) factorization |3l 0], which 
factorizes the production cross section into the production of a heavy quark pair in 
different non-relativistic states, multiplied by the transition for the pair to transmute 
into the observed quarkonium. The production of the heavy quark pair is effectively 
perturbative and is organized in powers of and v, the relative velocity of the heavy 
quark in the pair’s rest frame, while the corresponding transition rate is nonperturba- 
tive and is represented by a set of NRQCD long-distance matrix elements (LDMEs). 
If the factorization is correct to all orders in and powers in v, these LDMEs will be 
universal, i.e. process independent. Once the values of these LDMEs are extracted 
from one set of data, they can be used to explain and predict other data. 

However, the current NRQCD factorization formalism is far from perfect in de¬ 
scribing data on the heavy quarkonium production at high p^- The hxed-order 
NRQCD calculation suffers from large high-order corrections, due to the large en¬ 
hancement in the powers of and log(py/mQ)-type logarithms from high order 

contributions in powers of where tuq is the mass of heavy quark [5]. Although 
the most phenomenologically important contribution is assumed to be included in 
the NLO NRQCD calculation [6l [7] , the existence of large logarithm may potentially 
undermine the convergence of expansion at large pt 3> mQ. The lack of conver¬ 
gence may also be the reason for the long-standing heavy quarkonium polarization 
puzzle [Hinilin] and the recently-realized heavy-quark-spin-symmetry-violation puzzle 

[m la [13]. 

At large pt, the unstable perturbative series in the NRQCD factorization approach 
is due to the existence of another small parameter which is not included in 

the power counting of NRQCD effective theory. A satisfying framework for heavy 
quarkonium production at large pt must treat rn^Q/p^ in a systematic manner. 

Recently, a new QCD factorization formalism has been proposed to study heavy 
quarkonium production at large pt [a m [ig [ig nn ng. In this formalism, the 
production cross section is expanded in powers oil/px- It was argued to all orders in 
tts that the dominant leading-power (LP) terms, as well as the next-to-leading power 
(NLP) terms, can be systematically factorized into the perturbatively calculable hard 
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parts for producing a single parton (or a heavy quark pair at NLP), convoluted with 
corresponding single parton (or heavy quark pair) fragmentation functions (FFs) to 
the observed heavy quarkonium. Other than those suppressed by even higher powers 
in 1 /pt, all nonperturbative contributions to the cross sections are included in these 
FFs, whose scale dependence is determined by a closed set of evolution equations with 
perturbatively calculable evolution kernels. By solving the evolution equations, large 
perturbative log(p|./mQ)-type logarithms can be resummed to all orders in powers of 
as- Because of the systematic treatment of the powers of p'^/mq and log(py/mQ), 
the QCD factorization approach may converge faster in the perturbative expansion 
in powers of as than the NRQCD factorization. 


2 QCD factorization for heavy quarkonium 
production 

In the new QCD factorization approach, the production cross section of a heavy 
quarkonium H with momentum p at a large transverse momentum in the lab 
frame is expanded in a power series oi 1 /pt 0 [TH HH ITS] 
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where the hrst (second) term on the right-hand side is the LP (NLP) contribution, 
and the error of this factorization formula is suppressed by the power of 0{mQ/p'^) 
and higher. The LP (NLP) contribution is factorized into a sum of convolutions of the 
production of a single parton / (a heavy quark pair in state k) in the hard collision 
and the single-parton (heavy-quark-pair) fragmentation functions to all orders of as 
[5]. For the LP term, the parton / can be gluon, light quark, heavy quark or their 
anti-particles, while the QQ state k can be vector, axial vector or tensor, either color- 
singlet or color-octet for the NLP term. In Eq. ([I]), ^ is the longitudinal momentum 
fraction of the fragmenting parton(s), which is taken by the heavy quarkonium; Q 
and (2 are the relative longitudinal momentum of the QQ-pair in the amplitude and 
that in the complex conjugate of the amplitude, respectively, and can be different. 

In Eq. ([1]), the short-distance partonic hard parts da could be systematically calcu¬ 
lated in powers of as (we need to convolute with parton distribution functions (PDFs) 
if A and/or i? is a hadron). The FFs Df^niz] mg, fip) and Ci, C 2 ; rrig, pp) 
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are intrinsically nonperturbative, but process independent, universal functions. How¬ 
ever, their dependence on factorization scale /ij? is determined by a closed set of 
evolution equations [5], 
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^ ^[QQ{i^)]^\QQW)]^~i-' Ci> C2; Cl) C2)) 

where the evolution kernels 7 's and F's are process-independent and perturbatively 
calculable. The well-known DGLAP evolution kernels 7 /^/' are available to next- 
to-next-to-leading order in a^. The power-mixing evolution kernels 7 j_,.[qq(^/)] were 
calculated in Ref. |5], and the heavy quark pair evolution kernels have 

been recently calculated by two groups independently mu- If both K and k! are color 
singlet, the kernel r[gg(^)]_,.[gg(^/)] reduces to the well-known Efremov-Radyushkin- 
Brodsky-Lepage evolution kernel for exclusive processes [T9l|20]. By solving Eqs. ([ 2 ]) 
and ([ 2 ]), the perturbative logarithms log(p^/mg) are summed to all orders. 


3 Calculation of fragmentation functions with 
NRQCD factorization 

A set of input EEs at an initial scale /xq is required for solving the evolution equations. 

If the input scale is chosen as /xq ~ 2mg 3> AqcD) h is natural to use NRQCD 
factorization to further factorize these input EEs EZH 

^ ^ djgg(^jj_^jgg(^^j (Z, mg, PO) h'A)(Ci[gg(„)] (h-v))) 

[QQ{n)] 
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where /io and /ia are pQCD and NRQCD factorization scales, respectively. The 
summation of intermediate [QQ{n)] runs over all possible NRQCD states, which are 
labelled by spectroscopic notation j. The NRQCD LDMEs are ordered by 
powers of v, the relative velocity of the heavy quarkonium in its rest frame. In 
practice, the summation in Eqs. (01) and ([5]) can be effectively truncated and only 
hnite terms are kept on the right-hand side. 

With the NRQCD factorization, all the z, Q and C 2 dependences of input EEs 
can be calculated perturbatively and given by the short-distance coefficients d. The 
nonperturbative dynamics are absorbed in a few NRQCD LDMEs, which need to be 
determined by htting experimental data. We calculated these short-distance coeffi¬ 
cients to the hrst nontrivial order for all S'—wave and P—wave polarization-summed 
heavy quarkonium states. Conventional dimensional regularization was applied to 
regularize both UV and IR divergences. All derived short-distance coefficients d’s are 
IR-finite, and can be found in Refs. [211 122] • 

Since the polarization is an important observable for exploring the production 
mechanism of heavy quarkonium P El CD], it is critically important to study the 
production of polarized heavy quarkonia in the framework of the QCD factorization 
formalism, which requires a set of polarized input EEs. To make the calculation 
compatible with conventional dimensional regularization, in Ref. [23], we define the 
P—dimensional polarized NRQCD LDMEs by requiring that they preserve the same 
symmetry under rotation of angle vr about the pre-selected 2 ;-direction as their 4- 
dimensional counterparts. With the definition of P-dimensional NRQCD LDMEs, 
we calculate the short-distance coefficients in Eqs. (0|) and ([5|) to the hrst nontrivial 
order for all S'—wave and P—wave polarized heavy quarkonium states. The obtained 
d’s are all IR-hnite, and can be found in Ref. [23]. 

4 Comparison of LO QCD factorization with NLO 
NRQCD factorization for J/ijj production 

As shown in Fig. [H our LO factorized QCD contributions to the production rate [21] 
can almost reproduce the full NLO NRQCD calculation channel-by-channel for px > 
10 — 15 GeV. The comparison in Fig. [1] demonstrates that the very complicated and 
numerically-evaluated results of NLO NRQCD calculations can be reproduced by 
the simple and fully analytic LO calculation of the QCD factorization approach for 
Pt > 10 GeV, and clearly indicates that perturbative organization of the factorized 
power expansion is well suited to heavy quarkonium production at high px- It also 
shows the importance of the NLP contribution. Without it, as shown in Ref. [25], for 
example, the LP QGD factorization contribution can only reproduce NLO NRQGD 
results for the ^S'f^ and channels at large px, not for the and channels. 

To further illustrate the importance of NLP contributions, we plot the ratio of 
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Figure 1: Ratio of J/'0 production rate from LO QCD factorization over that of NLO 
NRQCD calculation for four leading NRQCD channels. See the text for details. 



Figure 2: Ratio of NLP contributions to total contribution in LO QCD for each 
channel, da means we have a special choice for PDFs. See text for details. 

the NLP contribution to the total LO QCD contribution in Fig. [2] for each channel. 
These ratios are independent of the values of the LDMEs. Figure [2] clearly shows that 
NLP contributions are negligible for the channel over the full pr range, and are 
small for the channel when > 20 GeV, beyond which it is below 10 percent. 
On the other hand, the NLP contributions are crucial for and channels, 
which remain very signihcant in the production rate, even if p^ approaches 100 GeV. 
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Since the FFs for a single active parton to fragment into a heavy qnark pair, 
calcnlated in NRQCD, vanish for both LO and NLO, the two-loop glnon FF derived 
in Refs. [26l [27] was nsed for the LP contribntion to the ^ channel in Fig. |2j 

In the above comparison with NLO NRQCD calcnlations, we did not include the 
evolution of FFs. A complete LO QCD calculation should include the evolution of FFs 
using the LO evolution kernels given in Ref. [S] and input FFs calculated in NRQCD 
factorization at NLO [TT] |22], and a set of updated NRQCD LDMEs by htting the 
data. From its consistency with the existing NLO NRQCD results, and the control 
through evolution of its higher order corrections, we expect the LO factorized QCD 
power expansion to do a better job interpreting existing data on heavy quarkonium 
production at collider energies, while a consistent NLO contribution is within our 
reach. Also, because the LP and channels, which produce predominantly 
transversely polarized heavy quarkonia, appear not to be dominant 01251 Eg, heavy 
quarkonium production at current collider energies is strongly influenced by the 
channel, and is more likely to be unpolarized. 
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